Spin density functional study on magnetism of potassium loaded Zeolite A 
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In order to clarify the mechanism of spin polarization in potassium- loaded zeolite A, we perform 
ah initio density-functional calculations. We find that (i) the system comprising only non-magnetic 
elements (Al, Si, O and K) can indeed exhibit ferromagnetism, (ii) while the host cage makes a 
confining quantum- well potential in which s- and p-like states are formed, the potassium-4s electrons 
accommodated in the p-states are responsible for the spin polarization, and (iii) the size of the 
magnetic moment sensitively depends on the atomic configuration of the potassium atoms. We 
show that the spin polarization can be described systematically in terms of the confining potential 
and the crystal field splitting of the p-states. 

PACS numbers: 73.22.-f, 75.75.+a, 82.75.Vx 



Searching for novel ferromagnets comprising only non- 
magnetic elements has been a fascinating challenge in 
condensed matter physics. Up to present, motivated by 
fundamental interests or potential technological impor- 
tance, various ferromagnets with non-magnetic elements 
have been synthesized. Among them, the alkali- metal- 
loaded zeolite is certainly a unique ferromagnet. Depend- 
ing on the crystal structure of the host cage and the num- 
ber/species of the guest cluster atoms, it exhibits not 
only ferromagnetism but also a rich variety of magnetic 
properties. In fact, one can envisage to design and con- 
trol the magnetic properties of this system by choosing 
appropriate combinations of the guests and hosts. ^ 

When we consider such materials design, the first step 
we should take is to clarify the mechanism of the spin 
polarization in this system. Indeed, it is of great interest 
to consider why clusters of alkali atoms confined in zeo- 
lite cages can be magnetic, even though the bulk alkali 
metals are usually non-magnetic. However, while almost 
20 years have passed since the seminal discovery of the 
spin-polarized ground state in zeolite A,^ ah initio stud- 
ies on magnetism in zeolitic materials have been quite 
limited so far. This is mainly because (i) the electronic 
structure is expected to be very complicated since the 
unitcell of the system is extremely huge [typically it has 
O(IOO) atoms and the lattice constant is ~ 30 A], and (ii) 
the atomic configuration has not been determined accu- 
rately in experiments due to the extreme complexity in 
the structure of the system. 

However, fortunately, the situation of potassium- 
loaded zeolite A (Ki2+nSii2Ali2048, hereafter we call it 
as K-LTA) is different; a detailed neutron powder diffrac- 
tion study has been performed^ and thus we exception- 
ally have a reliable reference for ah initio calculations. On 
top of that, recent density functional calculations based 
on local density approximation (LDA) have clarified that 
the low-energy electronic structure of this compound is 
quite simple.^ Namely, the system is regarded as a su- 
percrystal^ where guest potassium clusters act as a su- 
peratom with well-defined s- and p-like orbitals*^ Thus, 
we may expect that the low energy physics is expected 



to be described systematically in terms of a superatom. 

Experimentally, the spin polarization in K-LTA is ob- 
served when we introduce more than two potassium 
atoms (i.e., n > 2) and the Curie and Weiss temperature 
take their maximum at n ~ 4. If we follow the picture 
mentioned above, the superatom s state is fully occupied 
and three p states are partially filled for n ~ 4. While 
several scenarios for the spin polarization in K-LTA based 
on the assumption that the orbital degeneracy of the p 
states plays a crucial role have been proposed, ^^^'^^ an ah 
initio calculation based on spin density functional theory 
has yet to be done. The purpose of the present study is 
to clarify the mechanism of spin polarization in K-LTA 
from first principles. We examine in detail (i) whether 
a ferromagnetic ground state is really realized just from 
Si, Al, O and K, and (ii) if it does, for which condition 
the system has spin polarization. Hereafter, we focus on 
the case of n = 4, where the energy scale of magnetism 
is expected to be highest. 

Let us see the detail of the atomic configuration con- 
sidered in the present study. Figure [1] shows the atomic 
geometry determined by the neutron measurement with 
an assumption that the system has the symmetry of space 
group F23.^ Due to the Lowenstein rule (Si and Al should 
be arranged alternately), the unitcell contains two large 
(small) alminosilicate cages which we call a (/?) cages 
(upper left panel). While a simplified unitcell having one 
a and /3 cage (which violates the Lowenstein rule) is em- 
ployed in the previous LDA calculation^ here we adopt 
the original unitcell having ~200 atoms. Hereafter, we 
label the two a cages as ai and a2 (upper right panel). 
While the positions of Al, Si, and O were uniquely deter- 
mined by the measurement, the experiment found several 
possible positions for potassium atoms in the a cages 
(bottom panel). Following this observation, we denote 
the center of six-, four-, and eight-membered rings in the 
a cage as site I, II, and III, respectively. The site IV is 
the center of the a cage. In order to indicate which a 
cage the site I, III, and IV belong to, we label them with 
subscript such as site 1^^ or 1^2- Notice that the site II 
is on the border of two a cages. 
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FIG. 1: (Color online) Upper left: Overall profile of zeo- 
lite LTA, where O, Al, Si atoms are denoted by red-small, 
light-blue, and dark-blue spheres, respectively, a and p cages 
are marked by red-solid and blue-dashed lines, respectively. 
Upper right: Two kinds of a cages (ai and 0^2). Lower 
panel: Four kinds of sites occupied by potassiums (purple- 
large spheres). For the definition of the sites I-IV, see the 
text. 
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FIG. 2: (Color online) Four geometries considered in the 
present study, where dark-blue, light-green, and medium-red 
spheres stand for the potassiums occupying the site I, III, 
and IV, respectively. The difference in each geometry comes 
from the difference in the configuration of site-Ill potassiums; 
"two triangles" and "small triangle" (geometry I), "two tri- 
angles" and "large triangle" (geometry II), "hexagon" and 
"small triangle" (geometry III), and "hexagon" and "large 
triangle" (geometry IV) are accommodated in the ai and a2 
cages. We note that four site-I potassiums in the ai cage are 
overlapped with the remaining four site-I potassiums in this 



The experimental occupation numbers of site Iq,^, 1^2^ 
II, IIIc,,, IIIc,2, IVc,i, and IV«2 are 8.0, 8.0, 6.4, 5.9, 
3.5, 0.0, and 0.5, respectively. In the present calculation, 
for simplicity, we assume the occupation numbers as 8, 
8, 6, 6, 3, 0, and 1. Thus the site I, II and IVq;2 are 
fully occupied and the site IV^i is completely vacant. 
Note that the total number of the potassium atoms is 
32, and the number of site I, II, III in the unitcell is 16, 
6 and 24, respectively. While the atomic configuration 
of the site I, II and IV are uniquely determined, there 
are huge possibilities for the configuration of the site III; 
since there are 24 possible positions and we have nine 
potassium atoms in the site III (6 for the ai cage and 3 for 
the a2 cage), we have 12C6X 12(73=203280 possibilities. 

Thus in the present study, we focus on the following 
four geometries I-IV with trigonal symmetry shown in 
Fig. [21 In those geometries, we have two triangles and a 
hexagon in the ai cage, and a small triangle and a large 
triangle in the a2 cage. We note that, for the geometries 
I and III, the site-Ill and site-IV potassiums in the a2 
cage make a tetrahedron. 

Ab initio calculations based on local spin density ap- 
proximation (LSDA) were performed with Tokyo Ab ini- 
tio Program Package}'^ We used the GO A exchange- 
correlation functional^ plane- wave basis set, and the 
ultrasoft pseudopotentials^^ in the Kleinman-Bylander 
representation.^^ The pseudopotential of potassium was 
supplemented by the partial core correction. The en- 
ergy cutoff of plane waves representing wavefunctions is 



36 Ry and a 4x4x4 /c-point sampling is employed. 

We show in TABLE [H the total energies calculated for 
the geometries I-IV, together with local magnetic mo- 
ments in each a cage.^^ We see that the site-Ill configu- 
ration significantly affects the energetics and the value of 
the magnetic moment. Interestingly, for all the geome- 
tries, the moment develops in the a2 cage, which suggests 
that the valence electrons responsible for the magnetism 
reside mainly in the 0^2 cage. This is because the con- 
fining potential of the 0^2 cage is deeper than that of ai; 
as seen from Fig. [21 the site-I potassiums in the a2 cage 
exist at inner side than those in the ai cage, which makes 
the electrostatic potential in the 0^2 cage deeper. 

TABLE I: Total energy per unitcell of each geometry A^tot 
with reference to that of the geometry I and the local magnetic 
moments M^ in the ai cage {i — 1 01 2) (Ref. 17.V For the 
definition for the geometries I-IV, see Fig. [2l 





Geom I 


Geom II 


Geom III 


Geom IV 


A^tot(eV) 





1.02 


1.49 


0.12 


Ml (/is) 


-0.02 


-0.20 


0.06 


-0.16 


M2(/is) 


1.92 


0.33 


1.36 


1.00 



Another important point in TABLE [I] is that the size 
of the magnetic moment depends on the atomic config- 
uration of the potassium cluster. In order to clarify its 
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FIG. 3: (Color online) Calculated GGA band dispersions (up- 
per) and Bloch wavef unctions, ^a-^d, at the F point (lower) 
for the geometries I (left) and II (right). In the band disper- 
sion, red-solid and blue-dashed lines stand for majority and 
minority spin states, respectively. The energy zero is set to 
the Fermi level. In the wavefunctions, isosurfaces are drawn 
by the values of ±0.015 (a.u.). Potassiums occupying the site 
I, III, and IV are displayed by spheres of the same colors with 
Fig. El 



origin, we calculated the low-energy bands and their wave 
functions. In the upper left (right) panel in Fig. [3l we 
show the low-energy band structures near the Fermi level 
for the geometries I (II), corresponding to the second 
(third) column in TABLE [H While the band structure is 
quite different for these two geometries, the Bloch wave- 
functions at the r point (^a-^d) have common features; 
(i) they have large amplitude in the 0^2 cage due to the 
deeper cage potential in a2 mentioned above, and (ii) 
they widely spread in the a cage without localizing on 
any specific atoms and are regarded as superatom p-type 
wavefunctions.^ 

To make the situation clearer, in Fig.[4l we draw level 
diagrams of the low-energy states where the cage poten- 
tials, superatom p levels, and spin configurations are de- 
picted schematically. Here, the p states are represented 
in the local coordinates, where the z axis is taken to be 
perpendicular to the triangles in a cages (see Fig. [2]). In 
the geometry I (upper left), for both the a cages, the 
superatom px and Py levels are degenerated and higher 



than the Pz level. This can be understood as follows: In 
the ai cage, the superatom pz lobe extends in the di- 
rection of the center of the triangle and feels a positive 
crystal field, lowering the level of Pz. The same mecha- 
nisum works in the 0^2 cage, and the low- lying pz states 
in ai and a2 cages form a a bonding orbital. On the 
other hand, the Px and Py states in the a2 cage, with a 
non-bonding character, reside near the Fermi level. In 
the present case of n = 4, there are four superatom p 
electrons. Two electrons accommodated in the a bond- 
ing orbital and the remaining two electrons occupy the 
degenerated Px and Py orbitals. According to the Hund's 
rule, electron spins in the degenerated orbitals align to 
be parallel, thus generate the magnetic moment of ~ 2 
Mb. 

For the geometry II (upper right panel in Fig. [4j), in 
the a2 cage, the level of Pz is higher than px and Py. 
This can be also understood in terms of the crystal field 
formed by the site-Ill potassiums in the 0^2 cage. In this 
configuration, the Px and Py orbitals lie in the "large" 
triangle plane and thus these states are stabilized elec- 
trostatically, compared to the Pz state. Consequently, in 
the geometry II, the four electrons fully occupy the Px 
and Py states, so that the system has smaller magnetic 
moments. 

A similar argument can hold also in the cases of the 
geometries III and IV. In these geometries, the ai cage 
contains a hexagon instead of two triangles. The crystal 
field in the ai cage stabilizes the Px and Py states in the 
hexagon plane. In the geometry III (lower left panel in 
Fig.Hj), the two electrons occupy the low- lying pz orbital 
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FIG. 4: Level diagrams for electronic structures of the four 
geometries displayed in Fig. [2l where superatom p levels in 
the ai- and a2-cage potentials and the resulting spin con- 
figurations are drawn schematically. The p states are rep- 
resented in the local coordinates (for the definition, see the 
text). The difference between the cage potentials of ai and 
a2 comes from the difference in the configuration of the site-I 
potassiums (see Fig. [2]) and the p-level splitting is due to the 
site-III-potassium configurations. 
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in the 0L2 cage and the degenerated a orbitals made from 
the and Py states in the ol\ and cages accommo- 
date the remaining two electrons. Because the a orbital 
is less localized, the value of the moment is somewhat 
reduced (~ 1.5 /ie), compared to that of the geometry I. 
In the geometry IV (lower left panel), an extra state in 
the ^ cage appears near the Fermi level. This state has 
a large dispersion and small spin polarization. The net 
moment is generated from the three electrons occupying 
the degenerated 'p^ and Py orbitals, with the value of the 
moment ~ 1 /iB- 

Thus, we conclude that K-LTA can be ferromagnetic 
in LSDA and the mechanism of the spin polarization is 
consistently explained in view of the super atom picture; 
if the p levels formed by the superatom wavefunctions 
are degenerated and partially occupied, the system be- 
comes magnetic, due to the Hund's rule coupling in these 
states. This condition and the size of the moment are 
controlled by the atomic configuration of the guest potas- 
sium clusters, responsible for the cage-potential depth 
and the crystal field inducing the p- level shift and split, 
respectively. 

Experimentally, two kinds of possible magnetic struc- 
tures have been proposed for K-LTA. One is the ferri- 
magnetic model^ and the other is the spin-canted an- 
tiferromagnetic model. ^ In the former, the system con- 
sists of two sublattices with large (~ 2.8 /ie) and small 
(almost 0.0 /ie) magnetic moment, while, in the latter, 
every a cage has a moment of 1 \iB' These magnetic 
structures look like the results for the geometries I and 
IV, although the size of the moments is relatively smaller 
than the experimental values. It is interesting to note 
that the geometries I and IV have comparable energies 



(see TABLE III). 

In reality, the site III is randomly occupied.^ It is 
formidable or almost impossible to take into account the 
randomness of the site III explicitly in the first-principles 
calculation. However, the superatom picture has been 
shown to capture essential features of the low energy 
physics of K-LTA. Thus, a multi-orbital Hubbard-type 
model with random potentials would be served as a real- 
istic low-energy model describing magnetism of the sys- 
tem. Since all the parameters characterizing the Hubbard 
model can be derived with first-principle basis as was re- 
cently done in Ref. 18, we can discuss the magnetism in 
more realistic situations by solving this effective model. 
It will be an important future study. 

To summarize, by means of ah initio spin density func- 
tional calculation, we found that (i) K-LTA comprising 
only Al, Si, O and K can be ferromagnetic, (ii) the super- 
atom p states are responsible for the magnetism, and (iii) 
the magnetic property of K-LTA depends sensitively on 
the atomic configuration of the guest potassium cluster. 
The condition for the spin polarization was presented and 
understood systematically in terms of the superatom pic- 
ture. Our findings provide a firm basis and will be the 
first step for clarifying the low-energy physics in zeolitic 
materials. 
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